Universality in Glassy Low— Temperature Physics 
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We propose a microscopic translationally invariant glass model which exhibits two level tunneling 
systems with a broad range of asymmetries and barrier heights in its glassy phase. Their distribution 
is qualitatively different from what is commonly assumed in phenomenological models, in that 
symmetric tunneling systems are systematically suppressed. Still, the model exhibits the usual 
glassy low-temperature anomalies. Universality is due to the collective origin of the glassy potential 
energy landscape. We obtain a simple explanation also for the mysterious quantitative universality 
expressed in the unusually narrow universal glassy range of values for the internal friction plateau. 

PACS numbers: 05.20.-y,61.43.Fs,64.70.Pf,65.60.+a 



The physics of glassy systems at low temperatures dif- 
fers strikingly from that of their crystalline counterparts 
0. Differences are observed in thermal properties, in 
transport phenomena, and in dielectric and acoustic re- 
sponse, and are believed to be due to the existence of tun- 
neling excitations with a broad range of energy splittings 
and relaxation times, which are absent in crystals. The 
phenomenological standard tunneling model (STM) [2|, 
and its generalization, the soft potential model (SPM) [B|, 
which paraphrase this idea quantitatively were for many 
years considered to provide a satisfactory rationalization 
of glassy low-temperature physics. 

Questions into the origin of the universality of these 
phenomena Q| have not yet found clearcut answers; see, 
however, ||. In particular, the remarkable degree of 
quantitative universality (e.g., the narrow range of values 
for sound-attenuation in the 1 K regime) has remained 
a mystery so far. Also, relations between low- and high- 
temperature phenomena could not possibly emerge from 
phenomenological descriptions. Moreover, a number of 
recent experiments || Q || are difficult to reconcile with 
predictions of prevalent models. 

In order to better understand the successes and lim- 
itations of the phenomenological approaches, we have 
recently proposed microscopic models for glassy low- 
temperature physics |l0) using heuristics taken from 
spin-glass theory. The purpose of the present letter is 
to improve upon || [l0| with respect to both modeling 
and scope of analysis. We investigate a glass model de- 
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scribed by the Hamiltonian H = J2i 2m +^int({ui}) with 
a translationally invariant interaction energy of the form 
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The Ui designate deviations of particle coordinates from 
some reference positions, and (l|) may be thought of as 
arising from a Born von Karman expansion of the full 
interaction energy about these positions - taken to be 
stationary but not necessarily stable. As before || fTo[ . 
glassy properties are modeled by assuming the expan- 



sion coefficients at the harmonic level to be Gaussian 
random couplings with mean Jy = Jq/N and variance 
(SJij) 2 = J 2 /N. The non-random quartic potential is for 
stabilization, and so requires g > 0. Non-translationally 
invariant versions |{J or models with partial translational 
invariance JTo| have been introduced before. 

Our main results are the following, (i) The model ex- 
hibits a glassy phase at low temperatures, (ii) The po- 
tential energy landscape is represented self-consistently 
through an ensemble of effective single-site potentials 
which comprises both single-well potentials (SWPs) , and 
double- well potentials (DWPs), the latter with a broad 
spectrum of barrier heights and asymmetries, (iii) The 
distribution of parameters characterizing the single-site 
potentials depends on the collective state of the system, 
and differs qualitatively from what is assumed in the phe- 
nomenological models, (iv) The model exhibits typical 
low temperature anomalies, e.g., of the specific heat, and 
of acoustic or dielectric response, (v) Within the model 
we can explain, in particular the quantitative universality 
of glassy low-temperature anomalies expressed e.g. in the 
unusually narrow range of values for the internal friction 
plateau, (vi) We also see some unusual frequency depen- 
dence of internal friction curves similar to that recently 
observed in vitreous silica |^]. 

We now turn to the analysis. Due to the translationally 
invariant interaction, the system contains global trans- 
lations as zero modes, and the partition function must 
be evaluated orthogonal to these modes, i.e. with the 
constraint Ui — 0. Apart from this detail, the anal- 
ysis follows standard reasoning. The free energy is ob- 
tained by averaging an n-fold replicated partition sum, 
-j3f([3) = limAr^oo |tl _> (-Wn) -1 \n[Z^] av , and can be ex- 
pressed in terms of an Edwards- Anderson matrix of two- 
replica overlaps q a b- Although four-replica overlaps must 
also be introduced during the calculation, the final result 
does not depend on them. The replica free energy reads 
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with Z n a single-site partition function corresponding to 
the n-replica potential U n = |^ a (Jo + 12gq aa )u 2 — 
Eat 1a b u a u b -if(J2 a u 2 a ) 2 +gj: a u 4 a , that is, Z n = 
J Yla d u a exp{— f3U n }. The order parameters q a b must 
satisfy the fixed point equations q a b = {u a Ub) with angle 
brackets denoting a Gibbs average w.r.t. the potential 
U n . As usual, to perform the n — > limit, one starts from 
parameterizations of the q a t matrix based on assumptions 
concerning transformation properties of solutions under 
permutation of the replica. 

We have looked at the replica symmetric (RS) solution 
and a solution with one step of replica symmetry breaking 
(1RSB). Both in RS and 1RSB (and at all higher levels 
of Parisi's RSB scheme) the system is described by an 
ensemble of effective single-site potentials of the form 



U e s(u) = di u + d% u 2 + d^u 1 



(3) 



in which dn — g is the coupling constant of the quar- 
tic stabilizing interaction, and d\ and <i 2 are randomly 
varying parameters. We shall occasionally refer to d\ 
and d 2 as to an effective local field /i e ff and an effective 
harmonic coupling k e g via d-y = —h e s and d 2 = ^k e g. 
Their joint distribution is expressed in terms of the order 
parameters of the system which are in turn obtained in 
terms of Gibbs averages self-consistently evaluated over 
the U c s ensemble (|J). This ensemble of effective single- 
site potentials is a representation of the glassy potential 
energy landscape within a mean-field description. Glassy 
low-temperature anomalies follow from considering the 
effects of quantized excitations within this ensemble of 
local potentials. 

In RS one assumes q aa = qd and q a b = q for a ^ b. 
These must solve the RS fixed point equations 



<1 



(4) 



in which inner brackets denote a Gibbs average w.r.t. the 
effective single-site potential (||) with 



di = d? 



RS 



-Jy/qz 



d 2 = df = - ( Jo + I2g q d -J 2 C + J z) 



(5) 



and C = (3(qd — q), and outer brackets an average over 
the zero-mean, unit- variance Gaussians z and z in d\ and 
d 2 . The ensemble of single-site potentials thus comprises 
both SWPs, and DWPs, the latter with a broad spec- 
trum of barrier heights. A glassy state is signaled by 
q 7^ 0, thus a non-degenerate distribution of asymme- 
tries of single-site potentials. In the present model the 
transition is continuous and the critical condition is given 
by 1 = (J) c J) 2 ((u 2 ) 2 )_ which is just the de Almeida- 
Thouless condition for the occurrence of a RSB instabil- 
ity. For models with symmetric coupling distributions 
the transition temperature satisfies T c (J,g) = ^-T c (l, 1) 
with fc B T c (l,l) - 0.133 E* 111). In RS, d x and d 2 are 




FIG. 1: Central part of P(di\d 2 ), for J = 50 and g = 1. From 
top to bottom we have RS (independent of tfe) and 1RSB with 
d-i = 10, 5, 0, —5 and —10, respectively 



Gaussian and uncorrelated. This would be qualitatively 
in line with assumptions of the SPM, although there are 
quantitative differences. In contrast to the STM, corre- 
lations are predicted between asymmetries and tunneling 
matrix elements of tunneling systems in DWPs. 

The RS solution is unstable and thus strictly not 
acceptable throughout the glassy phase. For a non- 
translationally invariant model we have shown RSB 
effects to be small for the specific heat. Here, we look in 
greater detail also at the distribution of h c s and fc c ff. We 
shall find it to be qualitatively different from what has 
been obtained in the RS approximation, and thus also 
from what is assumed in the phenomenological models. 

In 1RSB, expected to exhibit the salient RSB effects, 
one assumes q aa — qd, q a b = qi for 1 < \a — b\ < m and 
qab = qo for \a — b\ > m. The fixed point equations are 
now much more complicated, and we shall not reproduce 
them here. Whereas d 2 has the same form as in RS, 
d msB _ i(j o + i2g qd - pC + J z), except that now 
C = f3(qd — qi), the local field is more complicated. It is 
formulated in terms of two variables z\ and zq 



d\ = d\ RSB = -Jsjq\ ~ Qo %i ~ Jy/^o z o 



(6) 



of which zq is a standard Gaussian, whereas Zi, while 
deriving from a Gaussian, becomes correlated with zq 
and z (thus d 2 ) in an intricate way through U c s- In the 
low-temperature limit, to which we restrict our attention 
here, the result is 



p(zi\z , z) 



exp{-z 2 /2 - D U eS (u)} 



2n J r>ziexp{-D U cS (u)} 



(7) 



Here D = (3m, which has a finite T — * 0-limit, u — 
u(z\, Zq,z) is the value of it which minimizes C/ C ff(u) for 
given values of zo, zi,and z, and T>z\ is a Gaussian mea- 
sure. Fig. 1 displays the resulting distribution of di 
conditioned on d 2 - 

The RSB result is qualitatively different from the RS 
result and from what is assumed in phenomenological 
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FIG. 2: Low T specific heat for J = 50 in RS (dashed) and 
1RSB (full). 



models, in that small d\ and thus symmetric tunneling 
systems are systematically suppressed. In fact, in the 
DWP region for negative d 2 and small d\, the form is 
P{d 1 ,d 2 ) ~ 7To(g? 2 )(1 + vri(d 2 )|di|), so P(d 1 ,d 2 ) is non- 
analytic at di = 0. Yet, due to the smallness of 771(^2)1 
RSB effects do not have significant influence on func- 
tional forms of thermodynamic and response functions 
at low temperatures. This was observed before for the 
specific heat in a non-translationally invariant model |)| , 
and is confirmed here. 

Glassy low-temperature anomalies follow from low en- 
ergy excitations in the ensemble of single-site problems 

2 

H c tf = f— + U c s(u). The energy of tunneling excita- 
tions in DWPs is determined by an asymmetry A and 
a tunneling matrix element Ao as E = \J A 2 + Aq, 
which arc in turn expressed through d\ and d 2 as 
A = E,d ly /2(\d 2 \ - 1) and A = £«v^|d 2 | 3/2 exp{|(l - 
M 2 | 3/2 )} H- Frequencies of higher (quasi-harmonic) ex- 
citations in SWPs or DWPs are given by the curvature 



of the potential in its minimum, U" s (u m i n j 
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By averaging over P(d±, d 2 ) one obtains tunneling and 
vibrational density of states and the specific heat as 
usual. Tunneling systems give rise to a slightly super- 
linear specific heat at low temperatures, C(T) ~ T 1+£ , 
with e decreasing with J (hence T c ), e ~ 0.1, 0.05 and 
0.01 for J = 25, 50, and 100 respectively, whereas the 
strongly peaked vibrational density of state is the origin 
of a Bose-peak in our model; see Fig. 2. RSB effects are 
small in C(T), the dominant effect being a reduction of 
the tunneling density of states as compared to RS. Ex- 
cept for the super-linearity at low T, which can be traced 
down to the translational invariance of the interactions, 
C(T) is qualitatively as in §, both in RS and 1RSB. 

It is interesting to compare with assumptions of the 
STM at the level of the distribution P(A, A) where A = 
ln(J5»/Ao). In the STM P STM (A,A) = P is assumed. 



We have 
P(A,A) 



P(d 1 ,d 2 



1 



3/(2|da|) 



(8) 



for d 2 < -(3/2) 2 / 3 , withdi =di(A,A) and d 2 = d 2 (A,X) 
instead. Unlike in the STM or the SPM (i) the result 
is not a constant, and (ii) symmetric tunneling systems 
are systematically suppressed. If one were to fit our re- 
sults to an STM parameterization one would associate 
P ~ P(0,A*) for some A* = 0(1) \d*\ = O(10). 
Interestingly, for reasonable energy scales (associated 
with tunneling of complexes like Si0 2 over atomic dis- 
tances, and giving between 1 and 5 K), we would esti- 
mate Pq = 10~ 7 . . . 10~ 6 Ar _1 /Atom for the system with 
J = 50, which is the right order of magnitude for typi- 
cal glasses. Also, using the fact that the low-T limit of 
the fixed-point equations entails qa, q oc J, and C cx J -1 
in RS, and similarly q<i,qi and qo oc J, C oc J -1 , and 
D oc J~ 2 in 1RSB, we can combine this with T c oc J 2 
and the expression for P(c?i,d 2 ) to predict the scaling 
Po ~ Tc 5/i for large T c both in RS and 1RSB, our sec- 
ond result relating low-T properties with a property of 
the glass-transition in our model. A detailed comparison 
with experiments is difficult, as a large change in glass 
transition temperature requires changing composition of 
the glas, and thus interactions in non-trivial ways. Yet 
the trend expressed by this result appears to be correct. 

Next, we discuss dynamics due to a coupling between 
local degrees of freedom u and the strain field e generated 
by Debye phonons, Psb — jue, where 7 is a deformation 

potential and e = ^Eg, s yj 2 ^7 ~ 6 L) the 

strain-field; the sum is over transversal and longitudinal 
modes and the tensorial nature of the strain-field is ne- 
glected. Quantities like internal friction are obtained by 
averaging the the dynamical susceptibility over the en- 
semble of single-site problems, 



772 X'L{u) ■ 



(9) 



At low T, and for huj <C k^T, the dominant contribution 
is the relaxational contribution of tunneling systems in 
DWPs within the ensemble of single-site potentials. For 
these a two-level approximation for H c g is appropriate 
and (|^) reduces to the well known expression known for 
two-level tunneling systems, when the internal coordinate 
u is approximated by a two state variable in terms of a 
Pauli matrix, u = uq<7 z . Fig. 3 shows the internal friction 
of the present model for J = 50 and driving frequencies 
as in m computed in 1RSB (except for a different global 
overall pre-factor, RS results are functionally hardly dis- 
tinguishable). Unlike in the STM and SPM, the domi- 
nant asymptotics is Or 1 - T 3+e (l + Cl T) for T < T* 

the 



, £ / 2 



andQ-^f^(^) (l + c 2 T)forT»T* 
influence of the c a = 0(1O -3 ) due to RSB effects being 
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supposition 7 ~ J expressing the fact that deformation 
potential and original interaction are of the same origin, 
and the elasticity theory scaling pv 2 ~ J for the sound 



T (a.u.) 



FIG. 3: Internal Friction as a function of T on a double log- 
arithmic scale for various driving frequencies, u) = 0.33, 1.03, 
2.52, 5.03 and 14.0 kHz (top to bottom). 



basically undiscernable. The exponent e is that intro- 
duced before to quantify the super-linearity of the low-T 
specific heat, and e ~ 0.05 for the J-value chosen. Note 
the following features: (i) there is a slight frequency de- 
pendence ui £ ' 2 of the plateau height, correlated with the 
low-T specific heat exponent (the extra T~ £ / 2 factor is 
undiscernable for the small range of temperatures con- 
sidered) Q ; (ii) the low T asymptotics is has by a tem- 
perature exponent slightly larger than 3. However, the 
crossover region from plateau to low-temperature asymp- 
totics is so large that effective exponents in accessible 
temperature ranges are still smaller than 3, the effect 
being stronger for lower frequencies. Both findings are 
well in line with recent experiments H , though observed 
effective exponents tend to show a somewhat stronger 
frequency dependence and be smaller than ours. 

Lastly, we turn to the universality issue. In uni- 
versality is explained within a renormalization group ap- 
proach as a collective effect due to interactions of quan- 
tized low-energy excitations. Here it is understood as 
a property of the interaction-generated glassy potential 
energy landscape, thus as a collective effect leading to a 
particular spectrum of quantized low-energy excitations. 
The mechanism is robust and may therefore justly be ex- 
pected to be insensitive to details. This holds in particu- 
lar for the tunneling density of states which is driven by 
a broad distribution of asymmetries and barrier heights. 
Higher order excitations depend to a larger extent on the 
shape of the stabilizing potential which is not collectively 
modified to the same extent, and so a stronger material 
dependence may be expected in this energy range; this 
fits well with the observed partial loss of universality in 
the Bose peak region. Within the present theory there is, 
beyond a rationalization of universality at the level of ex- 
ponents, also a rather simple explanation for the myste- 
rious quantitative universality of glassy low-temperature 
physics as expressed in the unusually narrow range of 
values for the internal friction plateau. It simply follows 
from combining the scaling Pq ~ J~ 5 / 2 with the plausible 



j-3/2 



T 



-3/4 



velocity, which together entail -j^- 
This weak parameter dependence due to cancellations 
implies, e.g., that the internal friction plateau is changed 
by only a factor 8 when interaction energies are changed 
such such as to increase T c from 100 to 1600 K! 

In summary, we have proposed and analyzed a transla- 
tionally invariant glass model. It exhibits typical glassy 
low-temperature anomalies of specific heat and acous- 
tic attenuation. The super-linearity of the low-T specific 
heat is linked with the non-degenerate barrier-height dis- 
tribution which in turn can be shown to follow solely from 
translational invariance (a feature that was absent in our 
original proposal |)|). Within the model we can correlate 
low- and high-temperature properties. Universality is un- 
derstood as consequence of collective effects, and we also 
have a simple explanation of the mysterious quantitative 
universality of internal friction data. 

This project profited a lot from a workshop on 
the low-temperature physics of glasses held at the 
MPIPKS-Dresden. Helpful discussions with C. Enss, U. 
Horstmann, H. Horner, S. Hunklinger, C. Picus and M. 
Thesen are also gratefully acknowledged. 
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